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Abstract. Wc study the dissipativc dynamics and the formation of entangled 
states in driven cascaded quantum networks, where multiple systems are coupled 
to a common unidirectional bath. Specifically, we identify the conditions 
under which emission and coherent reabsorption of radiation drives the whole 
network into a pure stationary state with non-trivial quantum correlations 
between the individual nodes. We illustrate this effect in more detail for the 
example of cascaded two-level systems, where we present an explicit preparation 
scheme that allows one to tune the whole network through "bright" and "dark" 
states associated with different multi-partite entanglement patterns. In a 
complementary setting consisting of cascaded non-linear cavities, we find that two 
cavity modes can be driven into a non-Gaussian entangled dark state. Potential 
realizations of such cascaded networks with optical and microwave photons are 
discussed. 
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1. Introduction 

The coupling of a quantum system to an environment is often associated with 
decoherence. This is exemphfied by the field of quantum computation, where the 
presence of additional reservoirs degrades the performance of quantum algorithms 
based on unitary operations executed on large many body systems [T]. On the other 
hand, in many quantum optical settings the environment is actually useful for realizing 
certain applications. Prominent examples are laser cooling or optical pumping, where 
quantum systems are prepared in highly pure states with the help of a reservoir 
[5], or continuous measurement schemes which enable the conditioned preparation 
of quantum states ^ . While quantum control schemes employing engineered unitary 
evolution are by now standard, the last decade has witnessed an increasing interest in 
alternative methods based on the concept of "quantum reservoir engineering", i.e., on 
controlling a quantum system by tailoring its coupling to an environment [4]. Efforts 
along these lines have lead, e.g., to a broad range of proposals for the dissipative 
preparation of entangled few body quantum states [51 |S1 [71 |H1 [HI [13 [H]- However, in 
recent years attention has in particular been devoted to the study of dissipative many 
body systems. In this context, it has been realized that quantum reservoir engineering 
allows one to dissipatively prepare interesting many body states [12l [HI [HI [151 [E] j 
perform universal quantum computation |14j . realize a dissipative quantum repeater 
|17j , or a dissipatively protected quantum memory [18j . Meanwhile, first experiments 
demonstrating the dissipative preparation of GHZ states in systems of trapped ions 
[T^ and EPR entangled states of two atomic ensembles [201 [H] have been reported. In 
these experiments the underlying principle has been to carefully design and implement 
a many-particle master equation, where a fully dissipative dynamics drives the system 
into a unique pure steady state representing the entangled state of interest. 

In this work, we study entanglement formation in driven few- and many-particle 
cascaded quantum networks as introduced by Gardiner and Carmichael |221 123j , where 
the unconventional coupling of multiple systems to a common unidirectional bath 
offers remarkable new opportunities for dissipative preparation of highly correlated 
states. As illustrated in figure [T] a cascaded quantum network consists of N systems 
coupled to a ID reservoir that has the unique feature that excitations can only 
propagate along a single direction, thereby driving successive systems in the network in 
a unidirectional way. While such a scenario is reminiscent of edge modes in quantum 
Hall systems [JH [25l [26l [27l [28], various artificial and more controlled realizations 
based on (integrated) non-reciprocal devices for optical t29j [SOj and microwave [211 122] 
photons are currently developed. Our goal below is to identify situations where 
cascaded quantum networks are driven by classical fields in such a way that they 
exhibit pure and entangled steady states. As shown in figure [ij such dark states 
emerge if the continuous stream of photons emitted by the first part ("A") of the 
network is coherently reabsorbed by the second part ("B"), such that no photons 
escape from the system and the output remains dark. The system hence acts as 
its own coherent quantum absorber while the constant stream of photons maintains 
entanglement between its two parts. 

We will first discuss the construction of coherent quantum absorbers on general 
grounds. Given the first part A of the network, we show how to choose the second 
part B such that the whole system evolves into a dark state. In doing so, we make 
use of the unidirectionality of the reservoir, which allows us to solve the first part 
independently of the second. These formal developments are then illustrated by two 
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Figure 1. A driven cascaded quantum network whicli is realized by a set of two- 
level systems coupled to a unidirectional bath. The two-level systems are driven 
by classical fields ^ Q and the continuously emitted radiation propagates along 
the waveguide and excites successive nodes. Under specific conditions, all photons 
emitted in a subsystem A are coherently reabsorbed in subsystem B. In this case 
the system relaxes into a dark state where no radiation escapes from the network, 
but a constant stream of photons running from A to B establishes entanglement 
between the two subsystems. 



settings in which the coherent quantum absorber scenario can be reaUzed. The first 
is a many body cascaded network, where each of the N nodes consists of a driven 
two-level system (TLS, "spin"). We show that this system exhibits a whole class of 
multi-partite entangled dark states, whose entanglement structure can be adjusted by 
tuning local parameters. As a second, complementary system we consider a network 
consisting of two non-linear cavities described by bosonic mode operators. By choosing 
appropriate laser drives and Kerr-type non-linearities one can ensure that this system 
evolves into a non-Gaussian dark state in which the two modes are entangled. In both 
examples, the coherent quantum absorber scenario thus leads to a dissipative state 
preparation scheme for non-trivial entangled states. In a more general context, these 
networks realize a novel type of non-equilibrium (many body) quantum system, which 
by changing the system parameters can be tuned between "dark" or "passive" phases 
(with no scattered photons emerging) and "bright" or "active" phases (with light 
scattered), while the nodes are driven into pure entangled or mixed states, respectively. 

The remainder of this paper is structured as follows. In section [2] we present 
the general model for an A^-node cascaded network and show how to construct the 
coherent quantum absorber subsystem B for some given subsystem A (cf. figure [I]) . 
Complementing these rather general developments, we discuss the two mentioned 
realizations of the coherent quantum absorber scenario in subsequent sections. Section 
|3] presents the driven cascaded spin-system for which we derive and discuss a multi- 
partite entangled class of dark states and also comment on the influence of various 
imperfections. Subsequently, section |4] discusses the setup based on cascaded cavities 
with Kerr-type non-linearity. In the latter case, we also obtain a purification of the 
steady state density matrix of the well-known dispersive optical bi-stability problem 
[33] as an interesting by-product. Implementations of the proposed cascaded networks 
are discussed in section [S] and concluding remarks can be found in section [6] 

2. Coherent quantum absorbers 

2.1. Cascaded quantum networks 

We consider the general setting of a cascaded quantum network as shown in figure [l] 
Here, N > 2 subsystems located at positions Xi are coupled to a ID continuum of right- 
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propagating bosonic modes b^, which represent, for example, photons in an optical or 
microwave waveguide. The whole network can be modeled by a Hamiltonian (h — 1) 



+ i/bath + ^ y" dc^ 5^ (46<,e"^^-/^ + H.c.) , (1) 



where -ffbath = / dw uj bl^b^^ is the free Hamiltonian of the bath modes and the integrals 
run over a broad bandwidth Aw around the characteristic system frequency loq. In this 
frequency range the bath modes are assumed to exhibit a linear dispersion relation with 
speed of light v. In equation ([T]) the Hi describe the dynamics of the individual systems 
and include classical driving fields (e.g. lasers or microwave fields). The system-bath 
coupling is determined by the "jump operators" Ci and coupling constants g^^, which we 
assume to be approximately constant over the frequency range Aoj. Implementations 
of an effective model of the form given in equation ([T]) can be achieved with atoms or 
solid state TLSs coupled to ID optical or microwave waveguides and will be discussed 
in more detail in section [5] below. 

The system-bath interaction in equation ([ij breaks time reversal symmetry and 
while photons can be emitted to the right, drive successive subsystems and eventually 
leave the network, the reverse processes cannot occur. To study the effects of this 
unconventional coupling, we assume that wg as well as the bandwidth Aw are large 
compared to the other relevant frequency scales and eliminate the bath modes in 
a Born-Markov approximation. This yields a generalized cascaded master equation 
(ME) for the reduced system density operator p [22l, f23| E] , 



.p-lY.[[c],c,p] + [pc],c,]y (2) 

i j>i 

Here the first part describes the uncoupled evolution of each subsystem Cip — 
—i[Hi,p] + 7l?[ci]p, where the Lindblad terms 'D[x\p = xpx^ — {x^x,p}/2 model 
dissipation due to emission of photons into the waveguide with a rate 7 — 2-Kg1^^. 
The unidirectionality of the bath is reflected by the last term in equation ([2|, which 
accounts for the possibility to reabsorb photons emitted at system i by all successive 
nodes located at Xj > Xi . The explicit Lindblad form of the ME ^ reads 

P = -i[i?casc, p] + l'D[c]p , (3) 

where i?casc — ^ ^2 Sj>i('']'^2 ^ '^I'^i) includes the non-local coherent part 

of the environment-mediated coupling, while the only decay channel with collective 
jump operator c = q is associated with a photon leaving the system to the rightjf] 
Note that equations ^ and ([3| are understood in a rotating frame in order to account 
for the explicit time-dependence of the driving fields. 



2.2. Coherent quantum absorbers 

In the following we are interested in steady state situations where every photon emitted 
within the system is perfectly reabsorbed by successive nodes in the network, such that 
there is no spontaneous emission via the waveguide output and the system relaxes to 
a pure steady state po = lV'o)(V'ol- To identify the general conditions for the existence 
of such states, we partition the network into two subsystems A and B as indicated 
in figure [l} with local Hamiltonians Ha and Hb, and jump operators ca and cb, 

X For N = 2 this ME stands in contrast to Ref. 1201 1211 . where the dynamics is purely dissipative 
with no coherent evolution. 
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respectively. Specifically, Ha = J2'iHi- i| Z]i>i(4'^» ~ '^I'^j) ^^'^ ^ where 
the primed sums run over the nodes of part A, and corresponding expressions hold 
for B. Then, in equation (j3|, iJcasc = Ha + Hb — i^(c^c^ — c^cs) and c = ca + cb, 
and the conditions for the existence of a pure stationary state are (see Ref. |13] and 



further comments in Appendix A I : 



(I) icA+CBMo)=0, (II) [ffcascPo] =0. (4) 

The first condition implies that the waveguide output is dark, i.e., (c^c) = 0, and the 
second one ensures stationarity. Within a quantum trajectory picture, condition (I) 
means that there are no stochastic quantum jumps, which would lead to a mixed state. 
In situations where [Ha^b,Po] = and ca,b|'0o) = for each subsystem separately, 
the network can simply be devided into two smaller parts which are then treated 
independently. In the following, we thus focus on situations where such a division is not 
possible and where the steady state possesses non-trivial correlations between A and B, 
as characterized for example by a non-vanishing C = (c^cb + cac'b) ~ 2Re{(c^) (cs)}. 
In view of (I) this third requirement can be expressed as 

(III) C = -2{{cicA)-\{cA)r)^0, (5) 

and directly connects the correlations between A and B with the amount of radiation 
emitted from the first subsystem. Note that for a pure steady state, C ^ also implies 
that A and B are entangled. Finally, we remark that under stationary conditions, a 
non- vanishing C implies a constant flow of energy from A to B, while the total scattered 
light vanishes. In the examples discussed below this 'coherent' absorption of energy 
in B can be understood as a destructive interference of the signals scattered by the 
two subsystems. 

The conditions (I) — (III) will not be satisfied in general. However, given a 
system A described by a Hamiltonian Ha and jump operator ca we can construct 
a perfect coherent absorber system B as follows. First, we point out that due to 
the unidirectional coupling the dynamics of A is unaffected by B, which can also be 
shown explicitly by tracing the ME ([3| over system B. In particular, the steady state 
of A is obtained by solving Cap'a — where Cap a = --'AHa^ Pa] + 1'D[ca\pa, and 
assuming a unique solution we write its spectral decomposition as p\ = '^i^Pk\k){k\. 
A pure state of the whole system is then given by |V'o) = Sfc ^/Pk\k)A ^ \k)B, where 
we assumed A and B to have the same Hilbert space dimension and defined |fc) = V\k) 



in terms of an arbitrary unitary V acting on Now, we demonstrate in Appendix 
|A]that conditions (I) and (II) can be satisfied by the choice 

(6) 

m\B, (7) 

where A„in = ("^|^^A,cff|^) and i?A,cff — Ha — i^c^c^i is the effective non-hcrmitian 
Hamiltonian associated with Ca, and we assumed a positive and non-degenerate 
spectrum {pk}- While equations ([6| and ^ define a general absorber system B, we 
find that for many systems of interest the stationary state satisfies ^/fhl{k\cA\n') = 
^Jpk{''^\cA\k) and ^Jjhi{k\HA.cs\'>'i') ~ \fVkk^\H a,cs\^ ■ In this case, the above relations 
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simplify to CB = -VcAVf and Hg = -VHAVf |] such that up to a unitary basis 
transformation the absorber system is just the negative counterpart of A. In particular, 
this situation applies to the two examples of cascaded spin systems and cascaded non- 
linear cavities, which we describe in more detail in the following sections. 



3. Cascaded spin networks 

Let us now be more specific and consider a set of N driven spins coupled to a 
unidirectional bosonic bath as shown in figure [T] In equation ([s]), the collective jump 
operator is now c = (t!_ and the cascaded Hamiltonian in the frame rotating at the 
frequency uJd of the external driving field reads 

Here the cr^ are the usual Pauli operators on site i, the fli are local Rabi frequencies, 
and the Si = uji — ujd are the detunings of the individual spin transition frequencies 
from the common classical driving frequency cod- In the following, the basis states of 
the spins are denoted by |e), \g), such that cr_ — \g){e\. The classical fields which are 
used to drive the spins can either be applied via additional local channels (assuming 
that the associated decay rate is much smaller than 7) or via a coherent field which is 
sent through the common waveguide. Note that by omitting the cascaded interaction 
in equation ([s]) we recover the familiar Dicke model for multiple two level atoms 
decaying via a common field mode, where even for O = a series of dark states can 
be identified from cjV'o) = 0. In contrast, in our cascaded setting these states are 
not stationary and non-trivial dark states can only emerge from an interplay between 
driving and cascaded coupling terms. 



3.1. Construction of dark states 

For N = 2 the dark state condition (I) restricts \4>o) to the subspace spanned by 
\gg) and the singlet \S) — {\eg) — \ge))/y/2. Condition (II) can then be satisfied for 
Qi = Q2 = ^ and any Si = —62 ^ S, for which we obtain the unique and pure steady 
state \^po) — \Ss), where 

\Ss) = ^=l={\gg)+a\S)), a ^ . (9) 

v/l + |ap 17-2(5 

The two spins thus realize a source and a matched absorber in the sense introduced 
in the previous section, and for the matrix representations of the various operators we 
can identify ca — o--, cb = —VcaV^ = o-_ and Hb = —VHaV\ using V = Gz- For 
strong driving, \a\ ^ 1, the state \Ss) approaches the singlet \S), where the mutual 
correlation |C| -> 1 is maximized. 

While for larger N a direct search for possible dark states is hindered by the 
exponential growth of the subspace defined by c\^q) = 0, we can use the state (|9| as a 
starting point and solve the cascaded system iteratively "from left to right" : Suppose 
that for Vti = Vl and Si = —S2 the first two spins have evolved into the dark state \Ssi) 
such that no more photons are emitted into the waveguide. Then, the following two 

§ We use this as a short-hand notation for an equivalent matrix representation of the two operators, 
{n\cB\m) = —{n\VcAV^\m) , etc. Here, the same set of states is used on both systems, which is 
permissible due to their equal Hilbert-space dimension. 
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Figure 2. (a) Circuit model for constructing the state |iS') in equation ( |13[ l for the 
example A' = {Sa, Si,, ... , —5a, —Si,, . . .). Each line represents a spin and is labeled 
by its detuning. A box connecting two lines denotes a unitary operation Ui{6i) 
and the corresponding exchange of detunings Si and <5i_|_i (see text and inset), (b) 
Scaling of the von-Neumann entropy of the first n spins in a network of A'^ = 12 
nodes for Q = 27. Upper curve: detuning profile as in (a) with \Si \ = -f/3. Lower 
curve: state corresponding to the detuning profile A' = (0, 7, —7, 7, . . . , —7, 0)/3. 



spins effectively see an empty waveguide and evolve into the dark state {Ss^), provided 
that ^3 = — ^4. By iterating this argument we see that for any detuning profile with 
621-1 — —S2i {i — 1,2, . . .) the steady state of the ME Q is given by 

= |55ji2 01553)34 , (10) 

which is shown explicitly in |Appendix B[ In particular, for the homogenous case 6i ~ 0, 
a strongly driven cascaded spin system relaxes into a chain of pairwise singlets. 

The dimer structure of the state 15") reflects the fact that radiation emitted from 
one node is immediately reabsorbed by the following one. We now consider situations 
where this reabsorption occurs by several of the following spins, leading to multi- 
partite entangled states. To this end, note that by starting from the state in equation 
(10 1 we can construct another dark state \S') = U\S'^) by any global unitary operation 
U with [U,c] — 0, while implementing the Hamiltonian H' = UH^ascU^ would ensure 
stationarity. However, for arbitrary unitary transformations H' will in general contain 
additional non-local terms, and wc must hence restrict ourselves to unitaries U under 
which Hcusc is form invariant. As an example, we write -ffcasc = ^casc(A), where 
A = (Si, 62, ■■■) is the detuning profile, and introduce the nearest-neighbor operations 



Ut{9i) = exp 



i— {a,+ai+i) 



oc exp 



1 — (7, 

2 



(11) 



where CTj; = (cr*,CT*,CT*). Then, by choosing tan(6'j 



Hr 



= (Si+i — Si)/ J we obtain 

:(A'), (12) 



with a new detuning profile A' = P^.i+iA, where Pi,i+i denotes the permutation 
of 5i and 5i+i (this is demonstrated in Appendix B). Thus, by starting from a set 



A'^ of alternating detunings as defined before equation (10 1, we can simply swap the 
detunings of nodes i and i + lto implement a new cascaded spin network with a unique 
stationary state \S') — Ui{9i)\S'^) . By repeating this argument, we obtain a different 
pure steady state for each permutation A' of A°. This class of states is given by 

|5') = ?7(A0^ A')|5"), (13) 
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Figure 3. Tuning a six-spin network througli briglit and dark states by adjusting 
local detunings. Tlie curves show purity and output intensity (c^c), wliere the 
detunings Si are interpolated linearly between the profiles given on the horizontal 
axis. Red arrows indicate the entanglement structure of the dark states. The solid 
lines show the ideal case, while for the dashed lines we included a finite on-site 
decay kq = 0.00257 (see section [3. 3| |. Parameters are f2 = 7, (5o = 7/5, 5i = 7, 
52 = 7/2, 53 = 0. 



where {/(A" — )■ A') is a product of nearest-neighbor operations Ui{9i), specified by 
the sequence of nearest-neighbor transpositions required for transforming A" into A'. 
A graphical representation of C/(A° — >■ A') in terms of a circuit model is shown in 
figure for a specific example. 

3.2. Discussion 

For large detuning differences \5i — 5i+i\ ^ 7 the unitary transformations given 



in equation (11) are SWAP operations ^1^ between neighboring sites. In this limit 
the states \S') remain approximately two-partite entangled, but with singlets shared 
between arbitrary nodes in the network. In contrast, for \5i ~ ~ 7 the 

Ui correspond to highly entangling VSWAP operations. Then, the entanglement 
structure can be much richer and in general the states \S') contain multi-partite 
entanglement between several or even all nodes. While in this case a full 
characterization is difficult, we point out that the Ui conserve total angular momentum 
such that the \S') approach multi-spin singlets in the strong driving limit. The 
amount of entanglement between subsystems now depends very much on the choice 
of the detuning profile, as can be seen from the two examples displayed in figure [2]jb), 
showing oscillating and linearly growing block entropy, respectively. More generally, 
the cascaded network can be driven into different types of states by simply adjusting 
local detunings. This is illustrated in figure [3j where an adiabatic variation of the 
detunings in a six-node network is used to prepare pure steady states with 2-, 4- 
and 6-partite entanglement, separated by "bright" (mixed state) phases where the 
conditions in equation (|4| are violated. 

For assessing the relaxation time that characterizes the approach of the system 
towards the presented steady states it is sufficient to examine the simple detuning 
profile leading to since the spectral properties of the Liouvillian are invariant 
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Figure 4. Influence of imperfections on the entanglement properties of the steady 
state, (a-d) show the steady-state concurrence in a two-node network for various 
imperfections, (a) On-site decays (see text), (b) Dephasing of the spins modeled 
by additional Lindblad terms Cp = Yj^T>[a\\p. (c) Waveguide losses, where 
ri is the fraction of photons that gets lost between the nodes (modeled by adding a 
factor — V to the last term in equation |2|). (d) Deviation from the asymmetric 
detuning condition modeled by a symmetric detuning offset 5i = 82 = e. In (a-c) 
we have chosen S\ = 82 = 0. (e) Influence of on-site decays on concurrences 
Ci,i+i of reduced two-spin density matrices for TV = 6 spins with 5i = 0. 

(f ) Influence of on-site decays on the steady state four-partite entanglement in a 
four-spin network, quantified by the measure of Ref. !36j , which is bounded by 0.5 
in this case. The detuning profile is given by A' = (0,-7/2,7/2,0). 



under the transformations (111. Numerical calculations for small systems suggest 
that the preparation time for these states scales efficiently with the number of nodes 
N . In particular, the uniqueness of \S^) implies that effects related to non-unique 
steady states found in related systems |35j are absent in our case. 



3.3. Imperfections 

Under realistic conditions various imperfections like onsite decays or losses in the 
waveguide can violate the exact dark state condition and the system then evolves to a 
mixed ("bright") steady state. This is exemplified by the dashed lines in figure |3]for 
the case of onsite decays, modelled by adding an additional term C p ^ KoJ^i '^[<^-]p 
to the ME ([3|. One clearly observes how the scattered intensity increases, while the 
purity drops as compared to the ideal case with kq = 0. To study the entanglement 
properties in such non-ideal situations, we show in figure |4f^a)-(d) the resulting steady 
state concurrence for various types of imperfections in a two-node network. In addition 
to onsite decays (panel (a)), we also consider intrinsic spin dephasing (panel (b)), 
waveguide losses (panel (c)), and small deviations from the ideal detuning profile 
(panel (d)). We see that different sources of imperfections lead to a qualitatively 
similar behavior and that in all cases the entanglement is quite robust and optimized 
for intermediate driving strengths fJ. 

For larger networks, the scattering of photons from the first nodes due to 
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Figure 5. (a) Two driven cascaded cavities with Kerr-type non-linearities. For 
simplicity, we assume the cavities to be driven through their back mirrors (red 
arrows), whereas the dominant decay channel 7 is provided by the front mirrors, 
(b) Level scheme for the construction of the pure steady state in the occupation 
number basis \nm) = |n} + |m}_ of the c± modes (A'' = n + m is the total photon 
number). Straight arrows denote Hamiltonian matrix elements, while wavy ones 
correspond to the decay of the symmetric mode. Blue dots indicate population 
in the dark state \19\. 



imperfections also affects successive spins, as shown in figure Qe) for a six-spin 
dimer chain \S^) in the presence of onsite losses. One clearly observes that the 
bi-partite entanglement in the dimers decreases with increasing number of previous 
nodes. To study the robustness of genuine multi-partite entanglement in the presence 
of imperfections, we employ the entanglement measure proposed in Ref. [36]. It can 
be evaluated in a straight-forward way and figure |4][f) displays the results for a four- 
partite entangled steady state in the presence of onsite decays. The behavior of this 
measure qualitatively agrees with the results for the concurrence in the two-spin case 
(cf. figure iga)). That is, we observe a tradeoff between the maximal achievable 
entanglement and the robustness of the state. Finally, figure [3] also shows that for a 
fixed N , different bi- and multi-partite entangled states are affected equally. 



4. Cascaded non-linear cavities 



The general concept of a coherent quantum absorber introduced in section |2] suggests 
that the formation of dark entangled states can exist also for systems other than spins. 
As a second non-trivial example where this can be shown explicitly, we now discuss a 
setting where the two cascaded systems A and B are represented by two Kerr non-linear 
cavities as depicted in figure [5ja) . The resulting distribution scheme for continuous 
variable entanglement is an alternative to other cascaded settings considered in this 
context [37]. We denote the two bosonic cavity modes by a and 6, and the dynamics 
of the system is governed by the ME ([3| with collective jump operator c = a -\- b and 
cascaded Hamiltonian 

H,^,,^HA+HB-(^{b^a-a^b). (14) 

Here, the Hamiltonian of the first cavity (system A, frequency (jJc,a) is 

Ha = ^a^a + Ka'^a^aa + iQ.{a^ - a) , (15) 

where we have already moved to a frame rotating at the frequency ojd of the external 
driving field, such that A = lOc.a — is the corresponding detuning of the cavity 
frequency and Q. the associated driving strength. Further, K denotes the strength of 
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the Kerr non-linearity. Motivated by our analysis of the cascaded spin system above, 
we assume that the Hamiltonian for the second cavity is given by 

Hb = -Ab^b- Kb^b^bb + in(b^ -b) . (16) 

Here we have chosen all constants to be identical to those used in Ha, such that the 
first two terms have opposite sign as compared to system A. As demonstrated below, 
this educated guess for Hb ensures that the cascaded ME exhibits a dark state. 



4.I. Steady state solution 

In order to show that the ME ^ indeed exhibits a pure steady state for the local 
Hamiltonians chosen above we exploit once more the conditions given in equation 
Q. From condition (I), i.e. (a + fe)|'0o) — Oj it is clear that the steady state should 
contain zero quanta in the symmetric mode{||] Therefore, it is convenient to change 
to symmetric and anti-symmetric modes c± = {a±b)/ \/2 and to write the dark state 
ansatz as l-^o) — |0)+ <8) |xo)-, with |xo) = J2n ^n\n), where |n) are Fock states in the 
occupation number basis. In order to exploit condition (II) we rewrite the cascaded 
Hamiltonian in terms of the symmetric and anti-symmetric modes, 

i7ca.c - iV2ncl + (a - i7/2 + K{N - 1)) + H.c. . (17) 



Here N — cl^C-^. + c_c^ is the total number of quanta, which is conserved by all terms 
except for those oc fl. Condition (II) is equivalent to i?casc|'0o) — M^o), and by 
projecting this equation onto jV-'o) we see that it can only be fulfilled for A = and 
hence 

n lA + 7/2 
-q;„_i , e = 7T7 , x= — . (18) 



n X + n — 1 iK iK 

This recursion is readily solved, such that the unique solution of conditions (I) and 
(II) is given by 

K^o) - |0)+ ® Ixo)- , with Ixo) = -J. £ -^^|n) . (19) 

N ^ Vn r(x + n) 

Tl— 

1 /2 

Here N = [0^2(2:, a;* ; 2|ep)] and 0^2 denotes the generalized hyper-geometric 
function |38j . Due to the simple tensor-product structure of jV'o) its only non- vanishing 
normally ordered moments are those of the c_-mode, i.e. 

ty.. (e*)"e'"r(x*)r(x) ^F^jx* + n,x + m-2\e\^) 

^^-'^'-^^ r{x*+n)T{x + m) ,F,{x* ,x;2\e\') ' 

and the moments of the original modes thus read 

Numerical studies suggest that there are no additional mixed steady states. 
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Figure 6. Characterization of the steady state \ipo) for K/'y = 0.01 (first row) and 
K/'y = 0.5 (second row). (a),(b) photon number in the first cavity normalized to 
resonant Unear response. (c),(d) hnear entropy 5iin of the first cavity as defined in 
equation \22\ . In all panels the dashed lines indicate the critical driving strength 
f!c, above which the semi-classical response becomes bi-stable (see, e.g., Ref. 1391 ). 



4-2. Entanglement properties of the steady state 

Despite the simple structure of the steady state lipo) in the c± representation, it is 
generally entangled when written in terms of the original modes a and b according to 



^ , X , X (-1)'" /(ri + TO)! , , 

Ti,m— 



where the a„ can be read off from equation (19). We characterize the entanglement 
with respect to this bipartition by the mixedness of the reduced state p'\ = 
trsll'f/'o) ("001} — BB^ of the first cavity, as measured by its linear entropy 

Snn = 1 - trMf} . (22) 

The results for two different strengths of the non-linearity are displayed in figure [6] 
For a better orientation, we show in panels (a) and (b) the photon number of the first 
cavity, normalized to the resonant response in the linear case, i.e. A = K = 0. Here, 
one clearly observes the well-known behavior of a single Kerr-non-linear cavity, which 

II As shown in Ref. 1131 a pure stationary state could also exist if \ipo) is an eigenstate of the jump 
operator, i.e. (a + fe)|V'o} = -^'IV'o) with A' £ C. However, in the present example this would imply 
that the symmetric mode is in a coherent state of amplitude A'. We do not expect this due to the 
non-linearity of the problem and it can indeed be shown that such an ansatz does not lead to a 
stationary solution unless A' = 0. 
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is characterized by a deformation of the response curve for increasing Q. Above a 
certain critical driving strength fi^ the classical response becomes bi-stable, while the 
quantum mechanical response curve exhibits a sharp step |33j . From panels (c) and (d) 
we see that the regions of pronounced non-linear response are those where the entropy 
is large, signaling a high degree of entanglement between the two cavities. In general, 
the features of the response are more washed out for larger non-linearities, where also 
the regions of high entropy are more extended. Note that the state |xo): a-nd hence 
also \ipo), is generally non-Gaussian, as can also be shown by evaluating appropriate 
measures [3D]. Finally, we remark that for i^T — ^ the steady state approaches a 
product of coherent states, i.e. jV'o) (X" | — with (3 = e/x, where |/3) 

denotes a coherent state of amplitude (3. We are thus able to recover the classical 
limit in which no entanglement persists between the cavities. 



4-3. The inverse coherent absorber problem 

Finally, let us take a different view on the analysis presented so far and emphasize 
once more that cascaded systems may be solved in a successive fashion "from left to 
right". This means that tracing the full cascaded ME ([s]) over the second subsystem 
B yields a closed equation for the first subsystem A, which in this case reads 

PA ~ I^aPa = — i[Aa^a + Ka} o) aa + iil(a^ — a), p^] + lT^\o\PA ■ (23) 

In section [2] we have assumed that the stationary mixed state solution of this equation 
including its spectral decomposition is known, which allowed us to explicitly construct 
a corresponding perfect absorber system B, such that the whole system is driven 
into a pure steady state. However, in many cases solving for the steady state p\ of 



equation ( 23 1 is a non-trivial problem by itself, which in the present case was first 
accomplished using phase space methods [33]. Note that by using an educated guess 
for the Hamiltonian Hb of the absorber system and then solving for the pure steady 
state iV'o) of the whole network, we have also indirectly solved the mixed state problem 
of a single cavity by computing p\ = trB{|V'o)(V'o|}- Its explicit matrix-elements in 
the Fock-state basis agree with the expressions found in the literature |3T], as do its 



moments ((a''')"(a)™) obtained as a special case of equation (20) 03]. The procedure 
presented above thus represents an elegant alternative way of obtaining the mixed 
steady state of the dispersive optical bi-stability problem. 

Avoiding to work with mixed states by using pure states on larger Hilbert spaces 
has found wide-spread use in the quantum information community, one of the reasons 
being that pure states can be manipulated more easily ]Jj . From this point of view, the 
steady state ji/'o) possesses additional relevance as a purification of the density matrix 
p\. This is an interesting result in itself, since obtaining such a purification generally 
requires knowledge of the spectral decomposition of p\ — a strong requirement even if 
p\ is known explicitly. In the above example we were able to circumvent this difficulty 
by directly constructing the purification I^/iq) as a steady state of the cascaded ME. 
Although so far this procedure is limited to systems where the corresponding absorber 
Hamiltonian Hb can be obtained by an educated guess^ it is intriguing to think about 
cascaded quantum systems also as an analytic tool for calculating stationary states of 
non-trivial open quantum systems. 
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Figure 7. Possible realizations of a single node of a cascaded spin network as 
shown in figure [l] (a) generic setup where the unidirectional coupling is achieved 
by a cavity connected to a circulator, (b) non-reciprocal superconducting circuit 
based on the proposal of Ref. |31| , (c) optomechanical transducer based on toroidal 
cavities [43lli4| . 



5. Implementations 

The two key ingredients for realizing cascaded networks as discussed in the previous 
sections are (i) the implementation of low-loss non-reciprocal devices for directional 
routing of photons, and (ii) achieving a coupling of single quantum systems to a ID 
waveguide that exceeds local decoherence channels. In the following, we discuss several 
implementations that fulfill these requirements, where a focus lies on optical and 
microwave setups. Requirement (i) is crucial for realizing the unidirectional coupling 
between the nodes and can in principle be fulfilled by standard circulators based on the 
Faraday effect (see, e.g., Ref. 22j for an optical implementation), but also by exploiting 
unidirectional edge modes in media with broken time-reversal symmetry [551 [2H 123 ■ 
However, in recent years there has been increasing interest in designing non-reciprocal 
on-chip devices which are integrable with, e.g., microwave circuitry j31l I32[ or nano- 
fabricated photonic components [221 (33 ■ Such elements would allow one to build up 
cascaded networks in a very controlled way in both the optical and the microwave 
domain, and thus constitute a promising way of implementing the scenarios discussed 
in this work. 

Concerning the second requirement, we first discuss the case of the cascaded spin 
network presented in section ^ Realizing a coupling of TLSs to a ID waveguide 
could, e.g., be achieved along the lines of the experiments reported in Refs. [151 
where superconducting qubits or atoms are coupled to transmission lines or hollow 
core optical fibers, respectively. However, one can also do without a direct coupling 
of the TLS to the waveguide by using the generic and flexible approach depicted in 
figure [tJ^ a). Here, the TLS is coupled to a cavity, whose output port is connected to 
a circulator or another non-reciprocal device. In the bad cavity limit k ^ g, where 
g is the TLS-cavity coupling, a series of these nodes results in the desired model 

with an effective TLS -waveguide decay rate 7 = 2g^/K. A particular realization 
of this scheme in the context of circuit cavity quantum electrodynamics |471 131j is 
shown in figure [7][b). Finally note that other systems like optomechanical transducers 
(figure [tJ^c)) have been proposed to realize a similar unidirectional coupling [ISISI]. In 
the case of cascaded cavities discussed in section |4l the most important requirement is 
the realization of strong Kerr-type non-linearities, which are of equal magnitude and 
opposite sign in the two cavities. Such a tunable non-linearity can be realized by a 
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dispersive coupling of the cavity field to a suitable four-level system, which has been 
analyzed both for optical gSl 132] as weU as superconducting microwave cavities [S^ . 

In summary, we find that the current experimental capabilities for realizing 
strong TLS-photon interactions in various systems, combined with the development 
of novel non-reciprocal devices in the optical and microwave domain, enable the 
implementation and design of various cascaded quantum networks. Here, the 
dissipative state preparation schemes described in this work could serve as an 
interesting application exploring the unconventional physical properties of such 
devices. 



6. Conclusions 



We have shown that photon emission and coherent reabsorption processes in cascaded 
quantum systems can lead to the formation of pure and highly entangled steady 
states. In the case of spin networks, this mechanism provides a tunable dissipative 
preparation scheme for a whole class of multi-partite entangled states, and instances 
of this scheme might serve as a basis for dissipative quantum communication protocols 
[T7] . For the case of two cascaded non-linear cavities we have identified a dissipative 
preparation scheme for non-Gaussian entangled states, and we have illustrated how 
cascaded systems can serve as a new analytic tool to evaluate the stationary states 
of driven-dissipative systems. More generally, our findings show that such driven 
cascaded networks realize a novel type of non-equilibrium quantum many-body system, 
which can be implemented with currently developed integrated optical systems or 
superconducting devices. 
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Appendix A. General construction of coherent quantum absorbers 



We present the arguments leading to the results quoted in section [Z2| Given a system 
A in terms of its Hamiltonian Hj^ and jump operator c^, we seek a suitable system 
B described by Hb and cb which perfectly reabsorbs the output field of A in steady 
state. In the following, we construct system B by requiring that the total cascaded 
system evolves into a pure steady state. The dynamics of the network is governed by 
the ME ^ with Hamihonian 

Hce.sc = Ha + Hb - i'^icAcl - c\cb) (A.l) 

and collective jump operator c ^ ca + cb, and assuming that c has no eigenvectors a 
pure state \iPq) is a stationary state of equation (|3| if and only if [T3] 

C\%)=0, [i?casc,|V'0)(V'0l] -0. (A.2) 



As discussed in section 2.2 the cascaded nature of the interaction allows us to solve for 



the reduced steady state of system A without knowing anything about B. Assuming 
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uniqueness, we introduce its spectral decomposition 

k 

and assume that the eigenvalues pk are positive and non-degenerate. A potential pure 
steady state |V'o) of the whole system can then be written as a purification of p^, i.e., 
\tpo) = J2k VPfcl^)^ ^ l^)-B Ul- Here, |fc) — V\k), and V denotes an arbitrary unitary 
and we assume subsystem B to have the same Hilbcrt space dimension as A. 

In the following, we write out the tensor-products more explicitly, such that the 
dark-state condition c\^po} — reads 

(cA(^l + l(^cs)|^o) =0, (A.4) 

with 

CA = '^{n\cA\m)\n) {m\ , cb = '^{n\cB\m)\n){m\ . 



To proceed, we plug the ansatz for j-^o) into equation (A.4), yielding 
J2 VPfe (^{n\cA\k)\n) (E) \k) + {n\cB\k)\k) (g) = , 

k,n 

and after relabeling the indices we obtain 

{VPk{n\cA\k) + y/p^Ck\cB\n)'^ \k) ® |n) = . (A.5) 

k,n 

Since the \k) form an orthonormal basis, this condition is equivalent to 

y/p^{k\cB\n) = -^/pk{n\cA\k) , (A. 6) 

which fixes the operator cb to be 

CB = - E W— (m|cyi|n)|n)(m| . (A. 7) 

To construct the Hamiltonian Hb for system B we exploit the second condition 

[-ffcasc, |V'o)(V'o|] — needed for a pure steady state, which is equivalent to i?casc|'0o) = 

X\tpo) for A G R. Note that cji/'o) = implies 

T t t T t t 

(cA «) - 8) cb)\^o) = ^2 (c^CA (8) 1 - 1 (8) c^cb)|V'o) , 

such that this condition reads 

(^HA,eS 1 + 1 ® i?lj,eff) I'Ao) = Mi^o) , (A.8) 

where we have introduced the effective non-hermitian Hamiltonians i?j,off — Hj — 
i^CjCj. Since a finite A would only lead to a global shift of Hb below, we can assume 
A = without loss of generality. We write the Hamiltonian as 

HB=J2{n\HB\m)\fi){m\, (A.9) 



and to determine the matrix elements in this expansion we start from equation (A 



and proceed as for the dark state condition equation (A.4) with the replacements 



ca — >■ Ha^cS Sind Cb — > Hg ^jj. As an intermediate result this yields 



{h\HB\m) = — {m\HA,cii\n) - i-(n|c^CB|m) . (A. 10) 

Pm ^ 
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To express the right-hand side of this equation fully in terms of operators on A, we 
employ the identity 



(n|c^CB|m) 



1 



PrnPri 



(to|capVaI") 



(A.ll) 



which can be derived with the help of equation (A. 6), and then make use of the 
stationarity of system A, 



(A.12) 



Equation (A. 10) then becomes 



1 



{n.\HB\m) = -- [ J^{m\HA,M + a/ — ("^l^lcffl^) . (^-13) 

which determines the Hamiltonian of system B. The expressions (A. 7 1 and ( A.13[ ) are 
the results quoted in section [2?2] of the main text. 



Appendix B. Dark states of cascaded spin networks 

We provide details regarding the cascaded spin network presented in section |3] The 
system is described by the ME ([s]) with collective jump operator c — ^^d i/casc 

defined in equation (|8|). 



Appendix B.l. Uniqueness of the steady state jiS") 



We show by explicit construction that the state \S^) given in equation (10 1 is the 
unique steady state of the network, provided that 821-1 = —821 and 57^ = for 
i — 1,2, . . .. To do so, we exploit the fact that the cascaded interaction allows for a 
successive construction of the steady state and start with the case N — 2. In this 
case, the steady state P2 is obtained by solving £^^^P2 = where — C12 is given 
by the block-wise Liouvillian 



-'),P 



7 



I?[c,,,+i]p(B.l) 



with I 



-o"L . We have already seen in the main text that a solution is given 



by P2 ~ l'^<5i)('^5i li and its uniqueness can be shown by calculating the characteristic 
polynomial of £12 and realizing that there is only one zero eigenvalue for 7 > 0. 

We continue with = 4 and write the ansatz for the steady state as p\ = 
\Ssi){Ssi \ ® p, where /i is a two-node density matrix. The four-node Liouvillian can 
be rewritten as 



c^^'^p = £12^ + £34^ - 7 {yli, ci2p\ + [p42> C34]) 



(B.2) 



and we note that Ci2\Ss-^){Ss-^\ = as well as ci2\Ss-^) — 0, such that the equation 
£'•^■'^4 = simplifies to £34^ = 0. However, this is just the two-node problem we 
have already solved and the unique solution for the four-node netwrok is thus given 
by pI = \Ssi){Ssi\ <8) liS^g) (553 1. By iterating this argument in blocks of two spins, we 
obtain the steady state jiS") given in equation (10). 
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Appendix B.2. Unitary form invariance of the master equation 

We briefly demonstrate that the cascaded Hamiltonian of the spin network is form- 
invariant under the unitary transformations (111 as stated in equation (12). In order 

as follows: 



to calculate UiHcascU^ , we write j 



i + 1 for brevity and rearrange 



-^casc ^ ^ 



.7 

5. 



(i 



2"^ 



Cl, 



nat 



k,l^i,j 

■IS 

k>l 



2^ + - 



where we have again used the piecewise jump operator c^^; = ci*: + . . . + ctL- Note 
that Ui(9) commutes with the first three lines and we thus focus on the last one. We 
abbreviate the operators appearing there hy A= {at, — cr^)/2 and B = a-'^a^_ — (J-'_a\. 
and observe that they transform into one another: 

(B.3) 
(B.4) 

Introducing the difference of detunings 5 — {6i — Sj)/2, a non-trivial form-invariance 
of the Hamiltonian -ffcasc under Ui{9) is realized if 

TUn\ s: A :7: 



U,{0)AUj{e) ^ A COS 20 + iB siii2e 
U,{0)BU}{0) = Bcos20 + iAsm20 



U^{0) 



5 A 



UK0) 



-5 A 



'2^ 



(B.5) 



That is, we require the detunings to swap and the cascaded part to remain invariant. 
It is easy to check that this requirement results in two equations, which are both 
solved for the choice tan 6* — —25/7 — ('^i+i ~ which was to be demonstrated. 
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